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Abstract. We consider the inverse problem to determine a smooth 
compact Riemannian manifold with boundary (M, g) from a re- 
striction Aj k of the Dirichlet-to-Neumann operator for the wave 
equation on the manifold. Here S and TZ are open sets in dM and 
the restriction As,n corresponds to the case where the Dirichlet 
data is supported on K + x S and the Neumann data is measured 
on M + x TZ. In the novel case where SDTZ = 0, we show that Ag n 
determines the manifold (AT, g) uniquely, assuming that the wave 
equation is exactly controllable from the set of sources S. More- 
over, we show that the exact controllability can be replaced by the 
Hassell-Tao condition for eigenvalues and eigenfunctions, that is, 

A 3 -<C||^||* a(5) , 3 = 1,2,..., 

where Xj are the Dirichlet eigenvalues and (<fo)^=i is an orthonor- 
mal basis of the corresponding eigenfunctions. 



1. Introduction 

Let (M, g) be a smooth, connected and compact Riemannian mani- 
fold with boundary dM. We consider the wave equation with Dirichlet 
data / e Cg°((0,oo) x DM), 

(1) {dl - A g )u(t, x) = 0, in (0,oo) x M, 

w|(o,oo)x9M = /, in (0, 00) x dM, 

u\t=o = d t u\ t=0 = 0, in M, 

and denote by w/ = u(t, x) the solution of ([!]). For open and nonempty 
sets S, 1Z C dM and T G (0, 00] we define the restricted Dirichlet-to- 
Neumann operator, 

feC™((0,T)xS). 
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Often we write A^ ^ = Al IgS 7l and = A^. When / is re- 

garded as a boundary source, the operator models boundary mea- 
surements for the wave equation with sources producing the waves on 
(0,T) x S and the waves being observed on (0,T) x 1Z. We consider 
the inverse boundary value problem to determine (M, g) from A^j § ^. 

A problem of this type is often called a complete boundary data prob- 
lem if S = 7Z = DM and a partial boundary data problem if S ^ DM 
or 1Z 7^ DM. A sub-class of the partial boundary data problems are 
the local data problems where S = 1Z 7^ DM. The inverse problems 
with local data and the analogous partial boundary data problems with 
S PI 7Z 7^ have been studied broadly and we will give below a brief 
review of this literature. On the contrary, problems with S H 1Z — 0, 
that is, problems with disjoint partial data, have remained open to 
large extent. We are aware of only two previous results: in jlQ] Rakesh 
proved that the coefficients of a wave equation on a one-dimensional 
interval are determined by boundary measurements with sources sup- 
ported on one end of the interval and the waves observed on the other 
end, and in [22] Imanuvilov, Uhlmann, and Yamamoto proved that 
the potential of a Schrodinger equation on a two-dimensional domain 
homeomorphic to a disc, where the boundary is partitioned into eight 
clockwise-ordered parts Fx, T 2 , . . . , r 8 , is determined by boundary mea- 
surements with sources supported on 5 = T 2 U T 6 and fields observed 
on 1Z = r 4 U r 8 . 

Inverse problems with partial boundary data are encountered in 
mathematical physics and in various applications. For example in med- 
ical imaging and in the geophysical imaging of the Earth, measurements 
can usually be done only a part of the boundary. Often it is not possi- 
ble to observe fields on the same area where sources are controlled. For 
example in oil exploration, explosives are used as sources and hence 
it is difficult to measure waves near the sources. Also, many inverse 
scattering problems, such as the transmission problems on a line, are 
equivalent to disjoint partial data problems. 

In this paper we consider the problem with S fl 7Z = and show that 
the inverse problem to determine (M, g) given A^ ^ has unique solution 
if the wave equation ([T]) is exactly controllable from S. We say that 
(M, g) is exactly controllable from S in time T > if the map 

(2) U : L 2 ((0, T ) x5)4 L 2 (M) x H~ 1 {M), 

U{f) = (u f (T ),d t u f (T )) 

is surjective. The condition by Bardos, Lebeau and Rauch gives a 
geometric characterization of exact controllability (21 [12] • In particular, 
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Figure 1. On left, the dogbone region with the Eu- 
clidiean metric satisfies the assumptions of Theorem [T] 
when S is the dashed blue part of the boundary and 1Z 
is the light red part of the boundary, see p2, Fig. 6]. On 
right, a detail of the dogbone region: for any point y on 
the light red arc, the point x is the closest point to y on 
the dashed blue part of the boundary We overcome the 
difficulties arising from non-convexity by using modified 



boundary distance functions in Section 4.4 



if M has a strictly convex boundary, then exact controllability is valid 
when every geodesic, continued by normal reflection on the boundary 
and having length T , intersects S. We refer to [5] for the precise 
formulation of the geometric condition in the case that DM is non- 
convex. 

For our purposes the exact controllability can be replaced by a spec- 
tral condition that is strictly weaker in terms of the size of S. Namely, 
let us denote the Dirichlet eigenvalues of — A g by 

< Ai < A 2 < A 3 •••->■ oo 

and the corresponding L 2 (M)-normalized eigenf unctions by <pj. That 
is, (0j,0fc)L2(M) = Sjk and 

-AgCpj = Xj(j)j on M, <f)j\dM = 0. 

We say that the manifold (M, g) satisfies the Hassell-Tao condition for 
eigenvalues and eigenfunctions with the set S C dM if there is Cq > 
such that for all orthonormal bases (0 J )°2 =1 of eigenfunctions 

(3) Xj < C 110^-11^2(5) , for all j = 1,2,.... 

We denote by d(x,y), x, y G M, the Riemannian distance function 
of (M,g). Our main result is the following. 
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Theorem 1. Let (Mi,g\) and (M2,g2) be C°° -smooth compact and 
connected Riemannian manifolds with boundary and let Sj C dMj and 
IZj C dMj be open non-empty sets with smooth boundaries for j = 1,2. 
Suppose that there are diffeomorphisms $ : «Si — > S2 and \1/ : TZ± — > IZ2 
and that (H) or (H ! ) holds, where 

(H) (Mj, gj), j = 1, 2, is exactly controllable from Sj in time T > 
and there is T > T + 2maXj =1 ^ maXj.^. d(x,lZj) such that 

^m, gi ,s u nJ = **(Am 2iS2i52 ^ 2 (^*/)) 

for allfeCZ°{(0,T)xSt). 
(IT) (Mj,gj), j = 1,2, satisfies the Hassell-Tao condition with 
the set Sj and 

AS lliriAfWl / = **(A^ AfWa (*»/)) 
for all f e C °°((0,oo) x S x ). 
Then (Mi, pi) and (M 2 , (72) ore isometric and there is such an isometry 
F : Mi -»- M 2 i/iai F| 5l = $ and = 

Notice that M mt is not assumed to be known a priori and a part of 
the proof is to construct M as a smooth manifold. The same is true 
for dM \ S~Un. 

While proving the above theorem, we develop a new geometric tech- 
nique that we call modified boundary distance functions. This tech- 
nique allows us to overcome difficulties arising from possible non-trans- 
versality between the geodesic flow and the boundary, see Figure 1. 
Such difficulties are present in various inverse problems, see e.g. [2S1 

sg. 

1.1. Previous results. The first results for inverse problems for the 
wave equation and the equivalent inverse problems the heat or the 
Schrodinger equations go back to the end of 50's when Krein stud- 
ied the one- dimensional inverse problem for an inhomogeneous string, 
u tt — c 2 {x)u xx = 0, see e.g. [30]. In his works, causality was transformed 
into analyticity of the Fourier transform of the solution. Later, in the 
late 60's, Blagovestchenskii P [TO] developed a method to solve one di- 
mensional inverse problems that was based on finite speed of wave prop- 
agation. In the late 80's Belishev [6] developed the boundary control 
method by combining the finite speed arguments with control theory 
to solve inverse problems for the wave equation in domains of the Eu- 
clidean space M ra . A turning point in the study of inverse problems for 
wave equation happened in 1995 when Tataru proved a Holmgren-type 
uniqueness theorem for wave equations with non-analytic coefficients 
[iol HH] . This and the earlier results on the boundary control method by 
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Belishev and Kurylev [8] solved the inverse problem for the wave equa- 
tion on a Riemannian manifold with complete data. For later develop- 
ment of the geometric boundary control method, see [21[7J[27l[28j|3T]. In 
all the above results measurements were assumed to be given either on 
the whole boundary or it was assumed that waves are observed on the 
same sets where sources are supported, that is, S = 1Z in our notation. 
For the wave equation, the problem where the closures of the source 
domain S and the observation domain 1Z do not intersect have been 
studied only in few papers, see [U] and the references therein, typically 
in the one-dimensional or radially symmetric cases. In [33] we studied 
the case when there are three disjoint domains £i,£2,£3 C dM and 
assumed that all three Dirichlet-to- Neumann maps Aj^j^, As 2 ,s 3 , and 
As 3 ,Si are known. 

The steady state solutions of the wave equation satisfy an elliptic 
equations and thus the inverse problems for elliptic equations can be in 
many cases considered as special cases of hyperbolic inverse problems 
with restricted data. A paradigm problem of this type is Calderon's 
inverse problem [13J that concerns the determination of the conductiv- 
ity from the elliptic Dirichlet-to- Neumann map. A smooth isotropic 
conductivity in a bounded domain of W 1 , n > 3, is determined by the 
elliptic Dirichlet-to-Neumann map as was shown in the seminal paper 
of Sylvester and Uhlmann [44J. In two dimensions the first unique 
identifiability result was proven in |38j for C 2 conductivities and in [3] 
the problem was solved for L°° conductivities. The corresponding in- 
verse problems for the elliptic Schrodinger equation has been solved in 
[TH H3] . For anisotropic smooth conductivity (corresponding to a gen- 
eral Riemannian metric) in Q C IR n , Calderon's inverse problem was 
solved in two dimensions in [13] using the isotropic result (38]. The 
needed regularity was reduced to L°° later in [1] . In the case of dimen- 
sion n > 3, Calderon's inverse problem is of geometrical nature and 
makes sense for general compact Riemannian manifolds with bound- 
ary, as was pointed out in [37]. This problem remains open, however, 
and we refer to [TS1 1201 EH] for partial results. 

The partial data problem for the isotropic elliptic equation in Q C 
]R n , n > 3, under certain geometrical restrictions on the sets S and 
7Z, has been solved by Kenig, Sjostrand and Uhlmann [29]. Also, the 
local data problem for anisotropic elliptic equations in Q C M 2 was 
recently solved by Imanuvilov, Uhlmann and Yamamoto [231 El] • For 
two-dimensional manifolds it was shown in [35] that the local boundary 
data for the Laplace-Beltrami equation with S = 7Z 7^ dM determines 
uniquely the manifold and the conformal class of the metric. Later 
the amount of needed measurements have been reduced in [121 I2D1 
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l21j 122] . For Riemannian surfaces the determination of the potential 
in a Riemann-Schrodinger equation with local data corresponding to 
S = 1Z 7^ dM has been solved in [16] . The above mentioned (25] is the 
only result with S fl 1Z = concerning elliptic equations that we are 
aware of. 

2. Outline of the arguments and notations 

We begin by showing in Section [3] that the condition (H) in Theorem 
[l] implies the condition (H'). That is, we show that the Hassell-Tao 
condition ^ is implied by the exact controllability ([2]), and that the 
operator is determined by the operator when 

(4) T > T + 2 max d(x, TV) 

and ([I]) is exactly controllable from S in time T . If we take S = dM 
and assume that M is non-trapping, then ^ was proved by Hassell and 
Tao [IT] without using exact controllability. In the erratum for |17j . 
it was noted that exact controllability yields ^ and this observation 
was attributed to Nicolas Burq. We give a short proof for this fact in 
Section [3] and show that ^ is strictly weaker than exact controllability 
in terms of the size of S. 

We will prove Theorem [T] in five steps that we will describe next. 
We denote SM := {£ G TM; \£\ g = 1}, that is, SM is the unit sphere 
bundle, and define d±SM := {£ G dSM; (f,Tf)p > 0}, where v is 
the interior unit normal vector on dM. We define the exit time for 
G SM \ d + SM, 

r M (x, := inf{s G (0, 00); 7(3; x, G dM}, 

where x,£) is the geodesic with the initial data 7(0) 
Moreover, we define the strip 

(5) J\f n := {{s, y) G (0, 00) x TZ; s< a n (y)}, 
a n (y) := max{s G (0,T M (y,u)]; d{j(s;y,u),TZ) 

and denote by the image of Af-n under the map 

(6) (s,y) j(s;y,v). 

The first step is to show that A^ determines (Mn,g). That is, we re- 
construct a piece of (M, g) in the boundary normal coordinates. Notice 
that a n (y) > 0, see e.g. [28] p. 50]. 

The second step is to show that A^ ^ and (Mn,g) determine the 
map 

A 5 ,b : / ^ M / |(o,oo)xfl, / e ^o°°(( > 00) x 5 )> 



= x, 7 (0)=e 
= *}, 
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where B C MJf is a small ball sufficiently far away from the topological 
boundary of M-r. The first two steps can be thought as bootstrapping 
steps that give us the interior data A^g. The remaining three steps em- 
ploy only interior data, and we will reconstruct the unknown manifold 
(M, g) by extending the known subset iteratively. 
We denote by x the center of B and define 

M B := (0, a B ) x S X M, 

a B := mm max{s G (0, r M (x, £)]; d(^y(s;x,^),x) = s}. 

We denote by Mb the image of Mb under the map 
(7) (s,0^l(s;x,0- 

The third step is to show that A^s determines (Ms,g). That is, we 
reconstruct a piece of (M, g) in the geodesic normal coordinates. Essen- 
tially the same geometric method can be used to reconstruct [Mn, g) 
and (M B ,g) from A^ ^ and A^, respectively. We will describe this 
method in Section HJ 

The fourth step is very similar with the second step. We show that 
A$^b and (Mb, g) determine A$ t B' for a sufficiently small ball B' C Mb- 
The fifth step is to iterate the third and the fourth step and to glue the 
local reconstructions together. We will give the details of the second, 
fourth and fifth steps in Section |5j 

3. Exact controllability and testing weak convergence 

of sequences of waves 

Let us recall that for open T C dM and / 6 Co°((0, oo) x T) we 
denote by w = u the solution of 

d 2 t u - A g u = 0, in (0, oo) x M, 

«|(o,oo)x9M = /, in (0, oo) x <9M, 

u\ t =o = d t u\ t =o = 0, in M. 

We extend the notation for open T C M int and / G Cq°((0, oo) x T) 
as the solution u* — u of 

d 2 t u - A g u = f, in (0, oo) x M, 

u\( ,oo)xdM = 0, in (0, oo) x dM, 

u\t= = 9 t u\ t=0 = 0, in M. 

Lemma 1 (Blagovescenskh's identity). Let T > and let F be open 
either in DM or in M mt . Then for 

^GL7 -((0,oo)x5), /GL7 co ((0,oo)xr). 
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we have 



(8) (u f (T),u^(T)) L 2 {M) = (f, (JA 5 , r - #A 5J ^J)^W)xr), 
where Rijj(t) := ip(T — t) and Jip(t) : = | fj* T 1 ip(s)ds. 



For a proof in the case T C <9M see e.g. [28, Lem. 4.15]. The case 
T C M mt is analogous. In the lemma, the Riemannian volume measures 
are used on (M,g) and on (dM,g). As we do not assume g\dM to be 
known a priori, the right-hand side of ^ for r = TZ is not trivially 
determined by the boundary measurement data A^. However, for 
our purposes it is sufficient to choose an arbitratry positive smooth 
measure dS on TZ. Then there is strictly positive /i 6 C°° (TZ) such that 



where dS is the Riemannian volume measure of (TZ,g). Moreover, 



where K := JK s ,n - RA S ^RJ. 

Let us assume for a moment that the exact controllability ^ holds 
and let (fj)f =1 C Cg°((0,oo) x TV). Then the functions u f ^(T ) tend 
weakly to zero in L 2 (M) as j — > oo if and only if the inner products 
g with / = fj and T = T tend to zero for all V e £ 2 ((0, T ) x 5). In 
this section we will describe a method to determine if (u^ 1 (T))'^ =l is 
weakly convergent under the weaker assumption Let us point out 
that the assumption ^ is needed only in this step. 

If TZ = S then Lemma jlj allows us to compute ||-u^ J (T) || i2 ( M ) an d 
we can easily determine if (m m ^'(T))°^ 1 tends to zero. In this case the 
assumption ^ is not needed since we can replace the exact controlla- 
bility by the approximate controllability, the latter of which holds for 
arbitrary (M,g) by Tataru's unique continuation [35], see e.g. Th. 
3.10]. The case TZ = S was solved originally in [27]. 

3.1. Exact controllability and the Hassell-Tao condition. We 

recall that the eigenvalues of the positive Laplace-Beltrami operator 
-A g with the domain H 2 (M) n H^(M) are denoted by < Ai < A 2 < 
A3 • • • — > 00 and the corresponding L 2 (M)-normalized eigenfunctions 
by <j>j. 

It is well-known, see e.g. p], that the exact controllability Q is 
equivalent with the continuous observability inequality, 



(9) 



dS = /j,dS, 



= (^(T),^(T)) L2(M) , 



(10) 



HOo, ™l)|| 



Hl(M)xL 2 {M) 



<C\\d u w 



L 2 ((0,T )x5) ' 
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where w is the solution of the wave equation 

0, in (0, T ) x M, 



(11) d^w-A g w 



W\(Q,T)xdM 



0. 



w\ t=0 = w , d t w\ t =o 



Wi 



in (0,T ) x dM, 
in M. 



Lemma 2. Suppose that (10) holds. Then there is C > such that 
(12) 



y/Xj + 1<C 



J 3\\L 2 (S) ' 



for all j = 1,2,... . 



by a larger time. Let Wq 
corresponding solution w in the eigenbasis 



Proof. Notice that if (j 10|) holds then it holds also when T is replaced 

bj and W\ = in (11). By writing the 

bj)JLi we see that 

(13) w{t, x) = COs(y/\jt)(j)j(x). 

Moreover, 

\\( w 0,Wi)f H x {M)xL 2 {M) = Uj\\ 2 H i {M ) = (#i,#i)i2(M) + 1 

= (Ag4>j, <f>j)L*(M) + 1 = Aj + 1. 



If T > 



1 



2VAi 
T 



then 



cos 



Xjt)dt 



T sin(2 v /A7T) T 
V " 2 



4VAi 



> 0. 



If also T > To, then we have by (13) and (10) that 

Xj + 1 = ||(w ,wi)||^i Wxi2(M) < C ||^w||^ ((0iT)xr) 



< C- 



\d u w 



f u; lli 2 ((o,T)xr) 



Jo 



cos z 



Ajt)c?t 



C / {d v (j>j) 2 dS. 



□ 



Example 1. Lei (M,g) be the Euclidean unit disc in M 2 and denote 

r s :={e 1 ™; 6e(0,s)}. 

Exact controllability from T s holds for s > 1/2 and does not hold for 
s < 1/2, see e.g. [5J. However, the condition (12) holds on S = T S for 
s > 1/4. 

Proof. An orthonormal basis of eigenfunctions can be chosen in the 
polar coordinates (r, 6) G (0, 1] x (0, 2it] as 



'mnl 



Cnm COs(m9) 



%n2 



Cn m sm(m8) 



where c nm is a normalization constant and z mn is the nth positive zero 
of the mth Bessel function J m . The corresponding eigenvalues are 
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Anmj = an d we have used the indices m = 0, 1, 2, . . . , n = 1, 2, . . . 
and j = 1,2. Notice that z mn ^ z m > n i if (m,n) 7^ (m',n'), see e.g. (13 
p. 484]. Thus any normalized eigenf unction corresponding to the 
eigenvalue z^ n can be written as = a(j) mn i + 60 mn 2 with a 2 + b 2 = 1. 
Moreover, for m > 0, 

H^ll^rr.) = cL^L^m^mn)) 2 / (a cos(m#) + bsm(m9)) 2 d6. 

Jo 

In particular, 

^ / m , l • / /m2^ -2(-l + (-l) m )a6 + mvr 

(acos(m^) + bsm(m8)Yd9 = — 

4m 

mix — 2 7r 1 /" 27r . . „. , , „.,, , A 

> > — = — / (acos(m9) + bsm(mO)) d9. 

~ 4m 12 12 L 



Thus the lower bound (12) holds for ri/4 with the constant C/12 where 
C is the corresponding constant for ri- □ 

3.2. Testing weak convergence of sequences of waves. 

Lemma 3. Let V C dM be open and nonempty and let 

T > max dfx, T). 

A sequence (uj)^ C L 2 (M) converges to zero weakly in L 2 (M) if and 
only if both (a) and (b) hold, where 

(a) For all sequences (^ m )™ =1 C C^°((0,T) x T) such that the se- 
quence (u^ m (T))~ =1 C L 2 (M) is bounded, there is C > satis- 
fying 

\(v h u^(T)) L2(M) \ < C for all I, m = 1,2, ... . 

(b) lim^ M (^(T)) L 2 (M) = for all $ G C o °°((0,T) x V). 

Proof. In the proof we omit writing L 2 (M) as a subscript. If * s 
weakly convergent to zero then (b) holds trivially and (a) holds since 
weakly convergent sequences are bounded. 

Let us assume (a) and (b). We will first show that (i>j)&, is bounded 
in L 2 (M). Let Vq G L 2 (M). Tataru's unique continuation implies 
approximate controllability, see e.g. [28, Th. 3.10]. Thus there is 
(ip m )™=i C C£°((0,T) x T) such that t^-(T) -»■ w in L 2 (M). Hence 
for all / 

|K«o)|= lim |(v,,^(T))| <C, 

m— >oo 

where C > is independent of I. Thus (^)^i i s weakly bounded, and 
hence bounded in the norm. 
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Let e > and fix m such that |i> — M^ m (T)|| < e. Then for large / 

I {vi , v ) | < sup 1 1 Vl \ I I k - v?" (T) 1 1 + (vi , (T) ) 
i 

< sup e + e. 
l 

As e > and Vq G L 2 (M) are arbitrary, converges weakly to 

zero. □ 

Let us reindex the orthonormal basis (0j)°^ =1 of Dirichlet eigenfunc- 
tions so that for j = 1, 2, . . . , the functions 

<f)jk, k = 1,2, . . . ,Kj 

span the space of eigenfunctions corresponding to the eigenvalue Xj. 
Here < Ai < A 2 < A 3 ■ • ■ — > oo and Kj is the multiplicity of Xj. Let 
us choose a positive smooth measure dS on S. Then there is a strictly 
positive function \x G C°°(S) such that ^ holds. As explained in 
pp. 5-6] the Fourier transform of the operator As_n with respect to 
the time variable is a meromorphic map, and its poles and residues 
determine the Dirichlet eigenvalues Xj and also the spaces 

Ej := span^-^-fcls; k = 1, 2, . . . , Kj} c C°°(S), 

for each j = 1, 2, . . . . 

Proposition 1. LetT > and (ip m )m=i C C^°((0, oo) x«S) and suppose 
that there is C > stzc/i i/iai 



Aj < C 1| ll^^ , fe = 1, •• • ,Kj, j = 1,2, . . . . 

TTien £/&e following are equivalent 

(i) T/ie sequence u^ m (T), m = 1, 2, . . . , is bounded in L 2 (M). 

(ii) For a// Ci > and {ej)f =1 C l7°°(5) satisfying 

(14) ej e Ej, WejW^g.^ < C iy ^X~, 

there is Gi > swc/j t/iat 

/(*) / ip m {t,x)ej(x)dS(x)dt) <C 2 



00 / /-T 



where Sj(t) := sin(y'Aj(T — t))/ y<Xj 



Proof. Let (i) hold and let ej, j = 1, 2, . . . , satisfy (14). As the choice 
of the orthonormal basis {<f)jk', k = 1, 2, . . . , Kj} in the jth eigenspace 
is unique only up to a rotation, we may assume after a rotation that 

fjzj = Cj d u (j)ji\s, j = l,2, 
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where Cj > are some constants. Moreover, the sequence (cj)^ =1 is 
bounded. Indeed, 



V A i f— 

Cj ~Q^~ — IM^Hl^S) = H^'Ul^S) - C II e 3 II L 2 (5;d5) — CCi^/\j. 

Denoting C := (CC Ci) 2 we have 

00 / r T r ~ \ 2 00 f r T f \ 2 



I 2 

Il 2 (AQ 



and (ii) holds. 

Let us now assume (ii) and let v G L 2 (M). We denote by Pj the 
orthogonal projection onto the j'th eigenspace. We may rotate again 
the basis {4>jk\ k — 1,2, ... , Kj} so that 

Pv 

(fiji = — — jj , for all j satisfying PjV 7^ 0. 

Then (v, cj)jk) = for all k > 2 and all j. 

We may choose := fi~ 1 d u <pji\s in (ii). Indeed, 

||a* ^f^'l II 1,3(5^5) < C \\dv<f>jl\\&(dM) - 

where the second inequality holds by [T7j. We have 

\{v,u^{T))\ 2 = if^,^) /% / ^dyfadSdtl 2 
■ =1 ./o ./s 



< 



3=1 3=1 KJo Js 



r-T 

^ ll w llL2(Af) su PXl ( / S J / ^mejdSdt 



and the sequence u^ m {T), m — 1, 2, . . . , is weakly bounded. □ 
3.3. Continuation of the data in time. 

Lemma 4. Lei S,1Z G dM be open and non-empty and suppose that 
7p is exactly controllable from S in time T . // T satisfies Q) t/ien 



A^^ determines A^? n . 
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Proof. Let 5 G (0, T ) satisfy 

T > T + 2 max g?(;e, 1Z) + 5, 

and let / G Cg°((0,T + 5) X 5). We choose ft, G C§°((0,25) x 5) and 
ft' eC™{{5,T + 5)x S) such that f = h + h'. As the coefficients of Q 
are time independent, we may translate in time and see that As^h'it), 
t G (5, T + 5), is determined by A^ K . 

We recall that A£ w : L 2 ((0,T)x5) -»■ H~\(0,T) xTZ) is continuous 
[52] . Let us suppose that ft G L 2 ((<5, T + <5) x 5) satisfies 

(15) As,nKt) = ^s,nKt), for t G (T + <J,T). 

Notice that such a function ft exists. Indeed, by exact controllability 
(§) there is ft G L 2 ((<5, T + 5) x 5) such that 

(u h (T + 5),d t u h (T + 5)) = (u l (T + 5), d t u l (T + 5)). 

As also h(t) = h(t) = for t > T + 5, we^have u h {t) = u h (t) for 
t > T + 5. We have shown that there is ft G L 2 ((5,T + 5) x S) 



satisfying (15). 



By ( 15 ) the Cauchy data of u h and u h coincide on (T + 5, T) x 1Z. 
Thus Tataru's unique continuation [35] implies that u (t + T + 5) = 
u h {t + T + 5) for t near (T - T - S)/2, see e.g. [281 Th. 3.10]. In 
particular, u h {t) = u h {t) for t > T — max x . eM d(x, 71). Analogously 
to the above case of Ag^h'it), we see that A s ^h(t) = A s ^h(t), t G 
(T,T + 5), is determined by A^. In particular, we have determined 

A s ,nf(t) = A s ,nh{t) + A s ,nh'(t), te{T,T + S). 

That is, we have shown that A^ n determines A^ 5 . By iterating the 
above argument we see that is determined. □ 

4. Local reconstruction of the manifold 

In this section we describe a method to reconstruct (M-ji,g) and 
(M B ,g) from A^^ and A^, respectively. We will first consider the 
local reconstruction method under the additional assumption that the 
functions a-ji and a B are known, and then show how these functions 
can be reconstructed from the data A^ and As,b, respectively. Let 
us recall that o~n{y) indicates the distance when the normal geodesic 
starting from y G 1Z hits to the boundary or to a point on the cut 
locus, see ([5]). The main difficulty when reconstructing an{y) is that 
the normal geodesic may intersect dM tangentially and this is hard 
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to detect from the data, see Figure 1. To deal with this difficulty we 
present a method that is based on "perturbing" the boundary dM. 

We will next give a series of lemmas that is common for the recon- 
struction method of a-ji and a B and that of {Mn,g) and (Mb,j). We 
start by introducing the modified distance function dh- Let T G M and 
h : T -+ R. We define 

40, y) ■■= d(x, y) - h(y), x G M, y g r, 

dh(a;,r) := inf d h (x,y), x e M, 
j/er 

where d is the Riemannian distance function of (M,g). Moreover, we 
define the modified domain of influence 

M(r, h) := {x G M; 40, r) < 0}, 

and denote for T > 

B(T, h- T) := {(*, y) G (0, T)xT; T — %) < t}. 

To simplify the notation, we define M(T, r) for a constant r G (0, oo) 
by M(r, h) where h(y) — r, y G T. Notice that if /i G C(T) then 

M(r, h) = {x G M; there is y G f such that d(x,y) < h(y)}, 

and our definition coincides with the definition of the domain of in- 
fluence in [39]. In particular, for T = {y} the set M(T,h) is the 
closed geodesic ball with radius h(y). In this case, we denote also 
M(y,h):=M(T,h). 

We will show first that As_n and determine certain relations 
between domains of influences, and then that these relations determine 
(M-ji,g) and (Mb,^). The latter step is purely geometric. 

4.1. From weakly convergent sequences of waves to relations 
between domains of influences. Tataru's unique continuation re- 
sult [15] implies that the wave equation Q is approximately control- 
lable, that is, we have the following lemma. 

Lemma 5. Let T > and suppose that V is open either in dM or in 
M mt and that h G C(T) satisfies h < T pointwise. In the case when 
T C M mt suppose, moreover, that h > pointwise. Then 

(16) {«/(T); feC?(B(r,h;T))} 

is dense in 



L 2 (M(T,h)) := {v G L 2 (M); supp(t;) C M(T,h)}. 
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Notice that the finite speed of propagation implies that u^(T) is 
supported in M(T,h) if / is supported in B(T,h;T). In the literature 
Lemma [5] is usually proved only in the case of a constant function h, 
see e.g. [28, Th. 3.10]. However, the case h G C(T) can be reduced 
to this case by approximating h with piecewise constant functions, see 
[391 Lemmas 4.2 and 4.3]. 

Lemma 6. Let T > 0, J G N, j — 1, . . . , J and let Tj be open either 
in dM or in M mt . Let hj G C{Tj) satisfy hj < T pointwise and, in the 
case Tj C M mt , also hj > 0. We define T := \J J j=1 Tj, 

(17) h{y) := ma.x{hj(y); j satisfies Tj 3 y} 

and denote Ui(f) '■— w(T). Let To be open either in dM or in M mt 
and let sq G (0, T]. Let fi G C°°(T U r ) be strictly positive. Then the 
following properties are equivalent: 

(i) M(r , S o)cU/=iM(r„^). 

(ii) For all f G C™(B(T , s Q ;T)) there is {f,)f =l C C?(B(T,h;T)) 
such that (Wi(/i(/o — fj)))^! tends to zero weakly in L 2 (M). 

Proof. Notice that / h- > jif is a bijection on C^°(B(T , s ; T)) and also 
on Co°(£>(r, h; T)). Thus we lose no generality by assuming that n = 1 
identically. We have 

M(T,h) = {x G M; there is y G T s.t. d(x,y) < h(y)} 
j 

= {jM(T u h 3 ). 

i=i 

The implication from (i) to (ii) follows from Lemma [5j We will now 
show that (ii) implies (i). We denote 

M :=M(T ,s ), M i: =M(r,/i), 
S :=B(T ,s ;T), S 1 := B(T,h;T). 

Let us assume that (i) does not hold and let x G Mo \ M\. As M\ is 
closed, there is a neighborhood U of x such that U fl Mi = 0. We will 
show next that U fl M™* is nonempty. 

If x G r , then points close to x are in M since So > 0. Let us now 
assume that x ^ T$. Then there is a path 7 from x to a point yo G Tq 
such that its length satisfies < £(7) < so. We may assume that 7 is 
a shortest path between x and yo and that it has unit speed pp. Then 
7(t) eUC\ M^* for small t > 0. 

We have shown that UdM^ is nonempty. Hence there is a nonempty 
open V C M such that V fl Mi = 0. By Lemma [5] there is a smooth 
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function fo supported in So such that j v u^°{T)dx ^ 0. However, by 
finite speed of propagation u$ {T)\ v = for any / supported in Si. 
Thus 

(u f °(T) - u f (T), l v ) = (« /o (T), l y ) ^ 0, 
for all / supported in Si and (ii) does not hold. □ 

Proposition 2. Let S, 1Z C dM and B C M mt be open, non-empty 
sets with smooth boundaries and suppose that the Hassell-Tao condition 
holds with the set S. Then As,n together with the smooth structure 
of S U 1Z determines the relation 

(18) 

{(?/ , yi,r, s, t) G ft 2 x (0, oo) 3 ; M(y , r) C M(K, s) U M(y u t)}. 

Moreover, A^^b together with the smooth structure of SUB determines 
the relation 

(19) 

{(l/o, 2/i, r, s, t) e OB 2 x (0, oo) 3 ; M{y , r) C M{B, s) U M(y u t)}. 

Proof. Notice that M(T,r) = M = M(T,T) for r > T and non-empty 
r C M since T > diam(M). The claim follows from Lemmas [lj[3] and [6j 
Proposition [T] and the following observation. Let yo,yi G M, so,si > 
and JgN. Let Tj be open either in dM or in M int and hj G C(Fj) 
for j = 1, . . . , J. Then the following properties are equivalent: 

(i) M{yo,ao)c\J J j=1 MF jt hj). 

(ii) For all e > there is a neighborhood r of 2/0 such that 

j 

M(r ,s ) C (jMO^ + e). 

If 2/0 G dM then we may take To C dM in (ii). Moreover, the following 
properties are equivalent: 

(i') M(2/ , so) C M(2/i, Sl ) U (j/=2 ^(Fj, hj). 
(ii') For all neighborhoods T 1 of 2/1 we have 

j 

M(2/ , s ) C M(r l5 Sl ) U |J M(Tj, hj). 

i=2 

If 2/1 G <9M then we may take T 1 C dM in (ii'). □ 
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4.2. Prom relations between domains of influences to distance 
functions. In this section we prove that the distances 

(20) d{j{s;y,is),z), {s, y) G M n , z G K, 



are determined by a-ji and the relation (18). Moreover, we prove an 



analogous result for the relation (19). To formulate the result, let us 
recall that we have defined for open T C dM and j/gT, 

a r (y) : = sup{s G (0, r M (y, u)); d(j(s;y,p),T) = s}. 

For open T C M int with smooth boundary and y G dT we define <Jr{y) 
by the same formula where v is now the interior unit normal vector of 
M \ r. We will show the following lemma. 

Lemma 7. Suppose that one of the following is satisfied. 

(a) T C dM is open and y G T. 

(b) T C M mt is open with smooth boundary and yo G dT. 

Letyi G M, t > and let < s < o"r(z/o)- Then the following properties 
are equivalent: 

(i) d(7(s;y ,v),yi) < t. 

(ii) For all e > there is 5 > such that 

M(y , s) C M(r, 8 -S)U M(y h t + e). 



In Section |4.4| we will reconstruct the function a^. To this end we 
will consider here also the modified distance functions dh- Let T C dM 
be open and let h G C 1 (r) satisfy 

(21) \grad 9M h(y)\ g < 1 for all t/GT. 

We define the modified normal vector, 

1 /2 

V(A) := (l - \gr&d dM h\l) v - grad SM /i, 

and the modified distance to a cut point 

a r (y; h) := sup{s G (0, r M (y, V(h))]; d h (y(s; y, V(h)), T) = s}. 

Notice that cr r (y) = a^iy; 0). 

To unify the notation we define V(h) and a-p(y,h) also for open 
T C M int with smooth boundary, y G dT and h = by the same 
formulas. That is V(h) = v and o"r(y; h) = crr(y) where v is now the 
interior unit normal vector of M \ T. 

Lemma 8. Let T C 6e open, h G C 1 (r), x G M mt and suppose 
that yo G T satisfies 

d h (x,y ) = d h (x,T). 
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Moreover, let 7 be a unit speed shortest path from yo to x. If h satisfies 
( lip then there is p > such that 7((0,p)) C M int and 7|[0,p] is the 
geodesic with initial velocity V(h). 

By |T] every shortest path is C 1 , whence we lose no generality with 
the assumption that 7 has unit speed. 

Proof. Let us denote t := d(x,yo). We prove the existence of p by a 
contradiction, so suppose that there is a strictly decreasing sequence 
(sj)^Li in (0, i) converging to zero such that 7(%) G <9M. Let us 
consider boundary normal coordinates (r, 2) G [0, 00) x DM in a neigh- 
borhood of yo. In these coordinates the metric tensor has the form 

n-l 

(22) g(r, z) = dr 2 + g°(r, z)dz 2 = dr 2 + 9%( r > z)dz j dz k . 

j,k=i 

We denote the boundary normal coordinates of 7(s) by (r(s), z(s)). 
Notice that for all e > there is 5 > such that for all s G [0, S] 

|zO)| s °(o,z( s )) < |z(s)l s °( 7 (s)) + e < l7( s )ls(7W) + e = l + e. 

Indeed, the first inequality follows from r(0) =0 and smoothness of 
g°, and the second one from (22). We denote yj := j(sj) and hi : = 
grad 9M /i. Then for small e > 0, large j and k > j, 



\h{yj) -h{y k )\ = I / d s h(z{s))ds\ < \\h'\\ c(T) / \z(s)\ g o {0z{s)) ds 

Js k Js k 

= \\h'\\c(T) i 1 + - S k) < Sj - s fc . 

By taking the limit — >■ 00 we have \h(yj) — h(y )\ < Sj. Hence for 
large j 

d h (x,yj) < Kl/\[sj,t]) - hfa) <t-sj- h(y ) + Sj 
= d h (x,y ) = d h (x,T), 

which is a contradiction as yj G T for large j. We have shown that there 
is p > such that 7((0, p)) C M mt . In particular 7|[0, p] coincides with 
a geodesic. The fact 7(0) = V(h) follows from a variation argument, 
see e.g. [36, p. 99] for a similar proof. □ 

Lemma 9. Let us suppose one of the following 

(a) T C dM is open, y G T and h G C(T) satisfies h(y) = 0. 

(b) T C M mt is open and has smooth boundary, y G dT and h = 
identically. 

Let t > 0. Then the following properties are equivalent 
(i) There is x G M such that d(x, y) = dh(x, T) = t. 
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(ii) For all s <t, M(y, t) <£ M(T, s + h). 

Moreover, if h £ C 1 (r) satisfies (21) and t < TM(y,V(h)), then x in 
(i) is unique and x = j(t; y, V(h)). 

Proof. It is clear that (i) implies (ii). Let us show (i) assuming (ii). 
We choose a sequence (xj,Sj) ! ^ 1 in M(y,t) x (0,t) such that Xj £ 
M(r, Sj + h) and Sj — > t. By considering a subsequence we may assume 
that X j T X • As M(y,t) is closed we have x £ M(y,t), whence 

dh(x,y) = d(x,y) < t. 

Moreover, Xj M(T,Sj + h) implies 

d h {x,y) > d h (x,T) = lim d h (xj,T) > lim sj = t. 

j—toc j— ¥00 

We have shown (i). 

Let us proceed to show the uniqueness. Let 7 be a unit speed shortest 
path from y to x. Then 7 coincides with y, V{h)) as long as it does 
not intersect dM. Indeed, this follows from Lemma [8] in the 
and from an analogous variation argument in the case (b). Finally, 
d(x, y) = t < TM^y, V(h)) implies that x = j(t; y, V{h)). □ 

Lemma 10. Let us suppose one of the following 

(a) T C dM is open, y e Y and h £ C 1 (r) satisfies h(y) = and 

(b) T C M mt is open and has smooth boundary, y £ dT and h = 
identically. 

Let r > 0. Then (i) implies (ii) where 

(i) t < a T (y; h) 

(ii) For allt £ (0,r] and s £ (0,t), M(y,t) <£ M(T,s + h). 
Moreover, if / ~f(-;y,V(h)) is transversal to dM then (ii) implies (i). 

Proof. If (i) holds then x = 7(t; V(h)) satisfies (i) of Lemma[9]for all 
t £ (0, r], whence (ii) holds. 

Let us now suppose that 7(-; y, V{h)) is transversal to dM. We will 
first prove that r > Tm{v, V{h)) together with (ii) yield a contradiction. 
Let us denote r := tm(v, V(h)) and let t £ (r, r). By Lemma [9] there 
is x £ M satisfying d(x,y) = dh(x,T) = t. Thus any shortest path 
7 from y to x coincides with 7(-; y, V(h)) on the interval [0, r]. By 
transversality 7(7") ^ TdM, whence 7 is not C l at r £ (0,t). This is a 
contradiction, since 7 : [0,t] — > M is a shortest path, whence it is C 1 , 
see pp. We have shown that (ii) implies r < tm (y, V(h)). 

Now we see that (ii) implies also (i) by applying Lemma [9] for all 
te(0,r}. □ 



20 



MATTI LASSAS AND LAURI OKSANEN 



Proof of Lemma^ We denote Xq := 7(3; yo,^)- By Lemmas 10 and [9] 
the point xq is the only point x G M satisfying d(x,yo) = d(x,T) = s. 
In particular, xq £ M(T, s — 5) for 5 > 0. 

If (ii) holds, then xo G M(y±, t+e) for all e > 0. Hence d(xo, y\) < t+e 
for all e > 0, and we have (i). 

Let us now assume (i) and let e > and x G M(y , s). If there does 
not exist 5 > such that x G M(T, s — 5), then d(x, T) > s — 5 for all 
5 > 0. Thus 

s > d(x,yo) > d(x, T) > s, 
and x = Xq G M(yi,t + e). We have shown (ii). □ 

4.3. Prom distance funtions to local reconstructions of the 
manifold. By Lemma [7] the distances, 

(23) d{j{s;y,v),z), {s, y) G M n , z G K, 



are determined by a-n and the relation (18). The considerations in [281 
Section 4.4.6] imply that the distances (23) determine (Mft,g) in the 
boundary normal coordinates ([6]). The reconstruction of (Ms,g) from 
the relation (19) is rather similar. However, we will describe it here for 
the sake of completeness. 

Let x G M and p > and denote B := M(x, p). Let us suppose that 
p is small enough so that B is contained in a normal neighborhood of 
x in M mt and consider the interior data A$ b- 

We let r G (0, p) and denote 

y r (0 := *y(r;x,£), ^ : dM(s,r). 

Then YJ. is a diffeomorphism and 7(s; = 7(s + r; x, £). Notice 

also that 



(24) r+ min a M (x,r)(2/) = o" 

y£dM(x,r) 



B 



Lemma [7] implies that Y r , a B and the relation (19) for B = M(x,r) 
determine the distances, 

tZ(7(s + r; x, f ), 7 (r; ac, £1)), £ , 6 e ^M, < s < ct b - r. 

Let us denote by dist the distance function d in the geodesic normal 
coordinates ([7]). We have shown that (B,g), a B and I^s,b determine 
the distances, 



(25) dist((s,rj),(r,Q), (s, V )eN b , (r,f) G (0,p) x ^M. 
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Lemma 11. Let (sq,T]o) G Mb and let us consider the differentiated 
distance function 

$(r,f) := d {StT ,- ) dist((s,r)), (r,f))l(»,»7}=(«o,t?o)> ( r '0 G (°>P) x 

Tnen there is ro > suc/i t/iat t/ie image of (0,r) x S^M under $ is 
open in S 1 ^ V0 )Mb for all < r < r . 

Proof. For y e M , n e and t > we denote exp y (tn) := j(t] y, if). 
As (soj^o) is i n a normal coordinate neighborhood of x = the same 
is true for (r, £) with r small enough. Thus 

$ tt (r,0 := grad (Si??) d^t((s,n),(r,0)|( s , f? )=(so,r ) o) = - Pex P(^, w )( r :0» 
where P is the projection, 



-Pf := p p, P : T^s ,rio)M'B — > S( SO!V0 )Mb- 



As exp, * x is a local diffeomorphism around and P is an open map, 

we see that $"((0, r) x S X M) is open in S( SO)Vo )Mb for small enough r. 
The claim follows by using the isomorphism S( So ^ Vo )Mb — > Sf Q V0 \Mb 
induced by the metric g. □ 



Lemma [IT] implies that the second order homogeneous polynomial 
g(s ,rj ) is determined on the space T? So ^Mb = ^ n by the distances 
(25). Thus g is determined also on TMb- To summarize, when we are 
given As t B, o~ B and (B, g) we can determine (Mb,j) in the geodesic 
normal coordinates ([7]). 

4.4. Reconstruction of distances to cut points using modified 
distance functions. In this section we show that for open V C dM 
the distance to a cut point o~r is determined by the relation 

{(y, t,h)eTx (0, oo) x C 1 ^); M{y, t) C M(r, fc)}. 

Moreover, we show that for a small ball B C M int the cut time <r B is 
determined by the relation 

{(y, t,s)eTx (0, oo) 2 ; M (y, t) C M(S, s)}. 

Lemma 12. Lei T C <9M be open and let y G T. TYien i/ie map 

or : r x c x (r) ->■ r 

is /ou>er semicontinuous at (yo;0). Moreover, if 0"r(yci; 0) < tm^/o?^) 
i/ien o"r is continuous at (yo;0). 
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Proof. We prove the semicontinuity by a contradiction, so suppose 
that there is a sequence ((yj,hj))JL 1 converging to (yo,0) such that 
lim inf j^oo or(yj] hj) < o~r(yo; 0). We denote ho — and 



a 



a r(yj] hj), Tj := T M (yj, V(hj)), for j > 0. 



As [0, cr ] is compact, we may consider a subsequence and assume that 
o~j —> o"oo- We let T G ((Too, o~o) and define 

Xj := 7(T; *i : = 4,0^, T). 

Notice that Xj is well defined for large j. Indeed, the exit time function 
t m is lower semicontinuous, see e.g. [TS], whence 



lim inf Tj ~> t$> o~q> T. 



We denote x := j(T] yo,iy). Then continuity properties of the mod- 
ified distances imply tj —> d(x ,T). Moreover, T > implies tj < T 
for large j, and T < o~ implies d(xo,T) = T. To summarize, we may 
consider a subsequence and assume that 

(26) tj < lim tj = d{x , T) = T < Tj. 

There is Zj G T such that dhAxj,Zj) = tj. By considering a subse- 
quence we may assume that Zj — > z^ G V. We see that z^ is a closest 
point to x in T since 

d(x ,z oo ) = lim d h Axj,Zj) = lim tj = T. 

However, T < cxo implies that ^ = yo, see e.g. [T4J pp. 144, 115]. By 
considering a subsequence we may assume that Zj G T since — > y Q G 



r. Lemma 8 and the inequality (26) imply that 



Xj = -y(tj;Zj,V(hj)). 

As T < o"o, the map (r,y) h-> 7(r;y, z/) is a local diffeomorphism at 
(T,y ) e (0,oo) x T, see e.g. pU p. 144, Th. III.2.2]. Moreover, the 
map 

a : C^T) x (0, oo) x T C^T) x M, 
a(h,r,y) := (h,j(r;y,V(h))) 

is a local diffeomorphism at (0, T, |/o) since its derivative is of the form 

Id 

A d( rt y)"/(r;y,v)\ r =T, y = yo 
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where A : C 1 (r) — > T X0 M is a continuous linear operator. In particular, 
there is a local inverse such that in a neighborhood of (0, T, y ) we 
have (3(h, j(r; y, V(h))) = (r, y). Hence for large j 

(tj, Zj) = f3( y h j ,j( y t j ; Zj, V{hj))) = f3(h j ,x j ) = p(hj, j(T; Vj , V(hj))) 



which is a contradiction with (26). We have shown that or is lower 
semicontinuous at (yo',0). 

Let us now suppose that a < r and show upper semicontinuity by 
a contradiction. To that end, we suppose that > <j . We let e > 
satisfy <r Q + e < min(<T 00 , r ) and denote Xj(e) = 7(cr + e;yj,V(hj)). 
Then for large j 

cr + e = d hj (xj(e),T). 
In particular, for small e > 

Co + e = lim 4 (^-(e^r) = d(7(a + e;y ,u),T), 

which is a contradiction with the definition of <Jq. □ 

Lemma 13. Let V C dM be open and let y G T. T/ien t/iere zs 
(Vjjhj)°? =1 C T x C 1 (r) converging to (yo,0) such that the geodesic 
s i—T- 7(5; yj, V(hj)) is traversal to dM, hj(yj) = and 

lim ariy^hj) = a r (y o ;0). 



In particular, 



liminf a r {y,h) = <T r (y ;0). 



Proof. By jTSJ Lem. 12] there is a sequence (yj,i]j)'jL 1 C d-SM 
converging to (yo,v) such that 7(-; %) is transversal to <9M and 
TM(Vj,Vj) converges to r M (y , v) as j — )■ oo. We may choose (hj) ! ^ =1 C 
C 1 (r) converging to zero such that 

Mfi) = and g rad aA/^(%) = V]\t*8M- 

In the case when o"r(?/o; 0) < r M (y , v), the claim follows immediately 
from Lemma 12 Let us consider the case cxr(2/o; 0) = r M (y , v). Then 

a r (yo; 0) = T M (yo, v) = liminf T M (yj, V(hj)) > liminf ortijj] hj) 
> o- r (y ;0). 

Moreover, by considering a subsequence we may assume that o~r(yj] hj) 
converges to lim inf j^^ o~r{yj', hj) as j — > oo. The second claim follows 
from the first claim and the lower semicontinuity of or at (yo; 0). □ 
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For open T C dM, y G T and /i G C(r) we define 

ar(y; h) '■— sup{t G (0, oo); t satisfies (ii) of Lemma [9 
Lemma 14. Let F C dM be open and y G T. T/ien 

liminf a T (y;h) = or(y Q ;0), 

where the liminf zs ta£;en over a// (y, h) G Tx C 1 (r) snc/i t/iat = 0. 



Proof. Lemmas 10 and 13 imply 



liminf (Tr(y,h)> liminf crpfy; h) — Privo] 0). 

(j/,ft)->(l/0,0) (2/,^)->(2/0,0) 



Let (hj)^t 1 be as in Lemma 13 Then by Lemma 



10 



liminf a r (y; h) < liminf ariyy, hj) = lim a v {yj\ hj) = o- r {y ; 0). 

□ 

Lemma 15. Lei x G M mi and let p > be small enough so that 
B := M(x,p) is contained in a normal neighborhood of x in M mt . For 
r > the following properties are equivalent: 

(i) r + p< a B . 

(ii) For all t G (0, r] , s G (0, t) and y G 95, M{y, t) £ M(B, s) . 



Proof. By Lemma 10 and (24) is enough to show that j(-;y,u) is 
transversal to dM if y G 95 satisfies r := r M (y, v) = mm z£ Q B r M (z, v). 
Let £ G S^M satisfy j(p; x, £) = y. Then 

is a closest point to y on dM. Thus 7(-; 2/ 0) v ) intersects <9M normally, 
in particular, it is transversal to dM. □ 

Summarizing, when we are given As^, we may use Proposition [2] to 
first determine a-ji by Lemma 14 and then to reconstruct the subset 
(Mn,g) by Lemma [7] and [281 Section 4.4.6]. Analogously, when B 
is as in Lemma [15] and we are given A$.b together with (B,g), we 
may use Proposition [2] to first determine a B by Lemma 15 and then to 
reconstruct the subset (M^,o) by Lemma [7] and Section 4.3 



5. Global reconstruction of the manifold 

In the previous section we have shown, under the assumptions of 
Proposition [2j that A s ^ determines (Mr., g). In fact, we have described 
a method to reconstruct (M#, g) in the boundary normal coordinates. 
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We will show next that Ag^ determines the interior data A$ t B for such 
a ball B C M-r that there is V C 1Z and to > satisfying 

BC M(T,t ) CM n . 

Let T > t . By the finite speed of propagation for the wave equa- 
tion Q, (Mft, g) determines u f (T) for all / G C °°((T - to,T) x T). 
Moreover, as the functions 

u'(T), feC™((T-t ,T)xT), 

are dense in L 2 (M(T,to)), we see using Lemma[l]that As^ determines 
u^(T)\b for if) G Co°((to, oo) x By varying T > t and noticing that 
the equation is invariant with respect to translation in time, we see 
that m^|(o,oo)x_b is determined for if) G Cq°((0, oo) x 5). That is, A S n 
determines A^b- 

Let x G M int and suppose that g is known in a neighborhood of 
.t. Then we can choose p > small enough so that B := M(x,p) 
is contained in a known normal neighborhood of x in M int . Let us 
suppose that A S ^ B is also known. We have seen that A^ ^ determines 
(Ms,g)- Let B' C M# be a ball. Then there is to > such that 

B' C M(B,t ) C M B . 

An argument analogous with the argument above shows that A^b de- 
termines As^B'. 

Let us now show how local reconstructions (Mb, g) and (Mg>,g) can 
be glued together. Let X\ G M B and x 2 G M B >. We have seen that the 
maps As t B and As,b' determine the maps As,B(xj,p), J = 1,2, for small 
p > 0. By Lemma ([I]) we can compute the inner products 

(u^(T),u^(T)) L2{M) 

= (u^(T),u^(T)) L2(M) - (^ 2 (T),^(T)) i2(M) , 

for G Cg°((0, oo) x B(xj,p)), and V G CS°((0, oo) x 5). The set 

67 p := {s G (0,oo); M{B(x 1 ,p),s) C M(E(x 2 , p), s)} 

is determined by Proposition [2j Moreover, x = y if and only if S p = 
(0, oo) for all p > 0. That is, we know how to identify points in 
Mb H Mb/. In particular, we can reconstruct the transition functions, 
whence we have constructed (Mb U Ms',g). 
Let us consider the collection 

M := {B; B C M int is a ball, 

As,n determines (Ms,g) and As,b}- 
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We have shown that Ai is nonempty and that 

if x G M B and B G M then M(x, p) G M for small p > 0. 
In particular, the open set 

U := (J M B C M int 
sex 

is nonempty. We will show next that it is closed in M int . Let x G M mt 
and let Xj G ?7, j G N, converge to x G M mt . Then there is a uniformly 
normal neighborhood W of x and j G N such that Xj G W 7 , see e.g. [351 
Lem. 5.12]. For small p > 0, B' := M(xj,p) G .M since Xj G M B for 
some 5 G M. As is uniformly normal, we have W C Mb>- Hence 
x G U and we have shown that U is closed in M mt . As M mt is assumed 
to be connected, we have U = M mt . 

We have shown that (M mt ,g) is determined by As^. The smooth 
Riemannian structure allows us to recover also the closure (M,g), see 
e.g. [2H1 P- 2116]. This proves Theorem [TJ 
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